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The  analytic  theory  for  ID  foil  currents  and  potentials  in  spiral  wound  cells  with  arbitrary  tab  config¬ 
urations  is  derived.  The  theory  is  extended  to  account  for  state  of  charge  gradients  along  the  length  of  the 
foils  revealing  how  tab  geometry  can  affect  the  slope  of  constant  current  voltage  curves.  Analytic  ex¬ 
pressions  for  cell  level  properties  are  derived  for  common  tab  configurations.  The  quality  of  many  tab 
configurations  is  compared.  A  straight  forward  method  is  described  for  the  common,  previously  unsolved 
problem  of  disentangling  the  foil  and  stack  contributions  of  a  measured  cell  impedance.  Many  other 
examples  of  practical  problems  and  solutions  are  presented.  The  analytic  theory  shows  how  to  directly 
calculate  the  foil  geometry  factor  introduced  in  a  previous  article  [1].  Simple  methods  are  presented  for 
adding  corrections  for  foil  effects  to  a  stack  centric  cell  model.  The  analytic  results  provide  a  lot  of  insight 
and  ease  of  use  that  numerical  solutions  cannot  easily  deliver. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Area  Specific  Impedance  (ASI)  is  a  useful  concept  for  charac¬ 
terizing  the  rate  or  power  capability  of  a  particular  cell  chemistry 
[2-4].  This  applies  for  both  AC  and  DC  impedance  techniques.  A 
long-standing  problem  with  using  ASI  analysis  for  larger  cell  sizes  is 
disentangling  the  impedance  caused  by  the  current  collection 
system  (tabs  and  foils)  from  the  stack  impedance.  Nelson  et  al.  [5] 
were  able  to  avoid  this  problem  by  building  cells  with  continuous 
current  collection  (current  flows  in  the  from  edge  along  the  full 
length  of  the  foils),  which  reduces  the  effective  impedance  of  the 
current  collection  system  to  near  zero.  However  it  is  not  easy  to 
mass  produce  cylindrical  cells  with  continuous  current  collection. 
Cells  with  a  limited  number  of  current  collector  tabs  placed  at 
strategic  locations  along  the  foil  length  represent  the  vast  majority 
of  commercial  cells  made  today.  It  would  be  very  useful  to  find  a 
method  for  subtracting  the  effective  foil  impedance  from  the  total 
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measured  cell  impedance  in  order  to  obtain  an  accurate  estimate  of 
the  stack  ASI.  We  will  show  how  this  can  be  done  for  any  tab 
configuration. 

To  our  knowledge  all  previous  attempts  to  model  the  effect  foil 
current  and  potentials  have  focused  on  using  numerical  methods 
[6-13].  In  this  work  we  will  show  that  to  a  very  good  approxima¬ 
tion  the  problem  of  ID  foils  with  randomly  positioned  tabs  can  be 
solved  analytically.  It  is  well  known  that  analytic  results  have  a  lot 
of  advantages.  Most  notably,  they  provide  much  deeper  under¬ 
standing  and  they  are  much  easier  to  use.  Others  [14-17]  have  also 
made  progress  with  analytical  solutions  for  discharge  capacity, 
electrolyte  salt  concentration  and  solid  phase  Li  concentration. 

The  field  equations  for  current  and  voltage  distribution  along 
the  length  of  the  foils  (x  direction)  have  been  described  previously 
[7].  All  foils  are  coated  on  2  sides,  which  means  that  the  standard 
commercial  spiral  wound  cell  is  actually  2  cells  in  parallel.  In  this 
work  we  use  the  symbol  2 is  for  the  total  stack  current  exiting  a  foil 
layer  and  the  current  entering  one  stack  layer  will  then  be  is.  We 
will  use  the  symbol  ps  to  represent  the  area  specific  resistivity 
[Q  cm2]  of  one  stack  layer,  so  the  resistivity  of  2  layers  in  parallel  is 
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V2Ps.  A  more  detailed  discussion  appears  in  Ref.  [7].  The  field 
equation  for  a  region  between  tabs  is: 

vjfo(x)  =  -2W(Ra+Rc)U(x)  (1) 

where  0s(x)  is  the  potential  across  the  stack,  is(x)  is  the  current 
density  [mA  cm-2]  running  through  one  layer  of  the  stack,  W  is  the 
active  width,  RA  =  pAlWAtA  is  the  linear  resistivity  [Q  cm-1]  of  the 
anode  foil,  pA  is  the  anode  foil  material  resistivity  [Q  cm],  tA  is  the 
anode  foil  thickness  [cm],  WA  is  the  anode  foil  width  [cm],  with 
similar  definitions  for  the  cathode.  Following  [7],  the  x  direction 
goes  along  the  length  of  the  spiral  wound  foils,  y  direction  goes 
through  the  stack,  z  direction  corresponds  to  the  height  of  the  cell 
or  the  width  of  the  electrode  foils.  Eq.  (1)  may  look  like  a  Poisson 
equation,  but  it  is  not  because  in  the  source  term  (RHS  of  Eq  (1 )),  is 
and  0s  are  intimately  related.  For  example  if  one  assumes  a  simple 
Ohmic  relationship  as  in 

0s (X)  =  Vo(x)  -psis(x)  (2) 

where  V0  is  the  open  circuit  potential  [volt]  across  the  stack  (which 
should  not  be  confused  with  the  open  circuit  potential  of  the  whole 
cell),  then  one  obtains 

(v2-«2)0s  =  -a2V0(x)  (3) 

where 


2.2.  Derivation  of  field  equations 

The  divergence  of  the  foil  current  must  account  for  current  flow 
through  the  stack  and  current  flow  from  the  tabs.  The  stack  portion 
is  smooth,  and  the  tabs  portion  can  be  approximated  with  Dirac 
delta  functions  [18],  <5(x).  This  is  approximate  because  the  tabs  and 
the  welds  between  the  tabs  and  foils  will  have  finite  width,  but  this 
width  is  very  narrow  in  comparison  to  the  length  of  the  foils.  Since 
the  Dirac  delta  function  is  not  commonly  used  in  the  electro¬ 
chemistry  literature  it  will  be  useful  to  review  a  few  properties  that 
will  be  required  for  understanding  the  developments  below.  One 
reason  for  developing  the  more  complicated  multisegment  solution 
in  Ref.  [7]  was  fear  of  the  delta  function.  All  of  the  properties  of  the 
delta  function  required  for  this  work  are  summarized  in  Table  1,  and 
the  delta  functions  and  its  first  two  integrals  are  plotted  in  Fig.  1. 
With  Table  1  under  our  belt,  fear  of  the  delta  function  becomes 
unwarranted. 

For  a  cell  with  an  arbitrary  number  of  tabs,  the  foil  current  di¬ 
vergences  are 

nA 

WXIA(X)  =  -2LWis(X)  +lJ2fj5(X-  Xj )  (5) 

.M 

VxIc(X)  =  2LWis(X)  -l£fj6(X-Xj)  (6) 

feC 


a2_2W(RA  +  Rc) 

Ps 

A  differential  equation  with  the  form  of  Eq.  (3)  is  commonly 
referred  to  as  a  modified  Flelmholtz  equation  or  a  screened  Poisson 
equation. 

Solving  for  the  foil  potentials  in  a  cell  design  with  multiple  tabs 
is  a  rather  complex  problem.  In  Ref.  [7]  we  handled  this  by  dividing 
the  foil  length  into  segments,  with  each  segment  bounded  by  a  tab 
or  the  end  of  the  foils.  Multiple  instances  of  (3)  were  numerically 
solved,  separately  for  each  segment.  The  boundary  conditions  for 
each  segment  are  related  to  the  external  boundary  conditions  for 
the  cell  through  a  system  of  linear  equations.  Solving  for  the 
segment  boundary  conditions  is  not  computationally  expensive, 
but  it  is  difficult  to  program  correctly.  Below  we  show  how  the 
general  problem  can  be  handled  with  just  one  segment,  which 
should  be  easier  to  understand  and  program. 

We  have  attempted  to  make  this  article  somewhat  compre¬ 
hensive  and  provide  a  number  of  expressions  for  the  various  cell 
properties  in  tables  and  appendices.  Casual  readers  can  jump  to  the 
results  section  to  see  how  tab  configurations  affect  cell  impedance, 
or  to  the  Worked  Examples  section  for  step  by  step  instructions  on 
how  to  use  this  technology  to  solve  real  world  problems. 

2.  Single  segment  theory 


where  I  is  the  externally  applied  current,  IA  and  Ic  are  the  currents 
along  the  anode  and  cathode  foils,  fi  is  the  current  fraction  flowing 
into  tab  j,  and  Xj  is  the  location  of  tab  j.  J2jeA  means  that  j  indexes 
only  over  the  anode  tabs  so  that  the  anode  and  cathode  each  have  a 
distinct  set  of  j s.  For  example  A  =  {1,2,3 ,...nA)  and 

C  =  {nA  +  1  ,nA  +  2,nA  +  3...nA  +  nc).  Since  these  summations  and 
their  integrals  will  appear  often  we  will  define 


k  =  'Em*-*) 

JSC 

^  =  Efj8(x~xj)- 

jeA 

X"  =  rAX\  +  rcX’i  =  ^2rjfjS{X-Xj) 


where 


r a  = 


rC  = 


Ra  +Rc 


jeA 


Rr 


Ra  +  Rc 


jeC 


(7) 

(8) 

(9) 


2.1  Reduced  units 

Before  proceeding  it  will  be  useful  to  use  dimensionless  units  in 
the  x  direction, 


where  L  is  the  total  active  length  in  the  x  direction.  We  will  also 
define  A  =  aL.  The  reader  is  reminded  that  (a2 /ax2)  =  V2  = 
(l/L2)(92/aX2)  =  (1/L2)V2. 


are  the  dimensionless  partial  foil  resistances  which  will  always 
satisfy  rA  +  rc=  1.  The  indefinite  integrals  E'  =  f  E"dX  and  double 
integrals  E  =  f  E'dX  are  then  easily  determined  from  the  relations 
in  Table  1.  The  fractional  currents  obey  the  normalization 
conditions 

Efj  =  Efj  =  1 

JeC  jeA 

Using  Ohms  law  we  can  convert  (5)  and  (6)  into  expressions  for 
foil  potentials 
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Table  1 

Properties  of  the  Dirac  delta  function  and  the  related  Heaviside  step  function. 


Property 

Formula 

Nomenclature 

Normalization 

J  5(x)dx  =  1 

Indefinite 

fd(x)dx  =  H(x) 

Heaviside  step 

Integral 

function 

Origin  shift 

J  d(x  -  x0)dx  =  H(x  -  x0) 

Shifted 
heaviside 
step  function 

Scaling 

Step  function 

*(«)  =  f 

H}(x)  =  JH(x-  Xj)dx  =  (x-Xj)H(x-Xj) 

Shifted  ramp 

indefinite 

integral 

function 

Ramp  function 

tirfW dX  -i(l-X,)2  0  <  Xj  <  1 

definite  integral 

For  an  Ohmic  system  we  can  use  (2)  to  eliminate  is  and  get  the 
modified  Helmholtz  equation  for  the  stack  potential. 

(v2  -  A2)*,  =  -A2V0  (x)  +/Rf2"  (11) 

where 


^  =  l2a2 


LWRf  _  Rf 

ips  “  Rs 


(12) 


We  have  introduced  the  total  stack  resistance  Rs  =  V2psILW  [Q]. 
Our  goal  is  to  find  an  analytic  solution  for  (11 )  and  use  the  solution 
to  calculate  cell  level  properties,  internal  current  distributions,  heat 
generation,  compare  the  effectiveness  of  different  tab  configura¬ 
tions  and  show  how  to  decouple  foil  and  stack  impedances. 


=  -LRAVXIA(X)  =  2  L2WRAis(X)-ILRAS" 
V2</>c  =  -LRcVxIc(X)  =  —2L2WRcis(X)  +  ILRCX'^. 


2.3.  General  solution  for  constant  A 

We  can  build  up  solutions  to  (11)  by  solving  for  the  Greens 
function  with  zero  slope  boundary  conditions 


These  can  then  be  subtracted  to  yield  the  final  field  equation  for 
the  stack  potential 

V^c  -  ^x4>a  =  vj&  =  L(Ra  +  Rc)[-  2LWU  +  IS"}  (10) 

where  L(RA  +  Rc)  is  the  characteristic  resistance  [Q]  of  the  anode 
and  cathode  foils.  Since  this  quantity  appears  almost  everywhere  in 
the  analysis  below  we  will  give  it  a  name  and  symbol 


Kf  =  L(Ra  +  Rc)- 


(v£-A2)0c  =  -8(X-Xj),  VX<PC  |x=0  =  Vx0G|x=i  =  0 

(13) 


which  has  the  solution 


0c  M) 


cosh  XXj  cosh  AX  +  cosh  AX^cosh  AX 
2A  sinh  A 


-Tsinhp|X-Xj|) 


(14) 


Rf  corresponds  to  the  total  impedance  of  the  anode  and  cathode 
foils  joined  end  to  end.  It  is  the  impedance  of  the  longest 
conceivable  foil  current  path  in  the  cell.  So  our  final  field  equation 
for  a  tabbed  spiral  cell  is 

v£</>s  =  Rf  [  -  2LWis  +  IS"] 


where 
X  =  1  -  X. 

Note  that  </>c(X,Xj)  is  perfectly  valid  when  tabs  are  at  the  ends  of 
foils  (Xj  =  0  or  Xj  =  1 )  which  means  that  all  tabs  can  be  treated  on  an 
equal  footing,  no  special  rules  are  required  for  terminal  tabs.  Fig.  2 
shows  a  plot  of  A2</>g(X,  1  /4)  as  a  function  of  X  and  A.  We  can  see  that 


-2.6 

-2.4 

-2.2 

-2.0 

-1.8 

-1.6 

-1.4 

-1.2 

-1.0 

-0.8 

-0.6 

-0.4 

-0.2 


Fig.  1.  Dirac  delta,  heaviside  step  and  ramp  functions. 


Fig.  2.  Greens  function,  A0|(X,Xj)  withX,-  =  1/4. 
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it  is  peaked  at  X  =  Xj  and  the  cusp  gets  sharper  for  larger  A.  This 
represents  the  response  of  the  stack  potential  for  an  infinitesimally 
narrow,  unit  current  source  at  Xj  =  1/4.  For  small  A  the  voltage 
response  is  spread  out,  and  for  large  A  it  is  localized. 

With  the  Greens  function  in  hand  we  can  build  up  the  solution 
for  arbitrary  source  functions  on  the  RHS  of  Eq.  (10)  or  (11).  For 
example  the  solution  to 

(v2-A2)</v  =  -X2V0(X).  (15) 

is  a  convolution  integral 

$V(X)  =  X2  J  $c(X,S)Vo(S)dS.  (16) 

The  complete  general  solution  for  constant  A  can  now  be  written 
as  follows 

ts(X)  =  MX)  -  iRf r(X)  (17) 

where 

m  =  J2r]fA(X’Xj)  (18) 


Vcell  =  AXc)  ~  MXa) 

=  rcMXC )  +  rAMXA)  ~  JRf(rcr(Xc)  +  rAr(XA) 

+  rArc[S(Xc)-X(XA)})  (25) 

which  looks  like  it  has  the  form  V0cv  -  /Rceii  with 

Reel!  =  Rf  (rcr(Xc)  +  vaT(Xa)  +  rArc[X(Xc)  -  s(XA)]).  (26) 

We  must  be  a  little  bit  careful  interpreting  rc</>v(^c)  +  as 

an  OCV.  This  quantity  represents  the  instantaneous  voltage  that  the 
cell  will  relax  to  when  the  current  is  switched  off.  After  running  the 
some  current  through  the  cell  the  state  of  charge  will  not  be  uni¬ 
form  along  the  length  of  the  foils.  In  this  case  when  the  external 
current  is  turned  off,  there  will  still  remain  some  internal  mixing 
currents  flowing.  This  will  cause  the  external  cell  voltage  to  further 
relax  on  a  longer  time  scale. 

Up  until  this  point  we  have  expressed  all  quantities  in  terms  of 
the  characteristic  foil  resistance  Rf.  What  happened  to  stack 
impedance  in  Eq.  (26)?  Well,  it  is  hiding  in  the  T  summations,  Eq. 
(18).  It  is  also  useful  to  note  that  because  of  Eq.  (12),  for  any  of  the 
derived  expressions  in  this  article  we  can  substitute  RSA2  for  Rf.  We 
prefer  to  use  Rf  because  it  can  easily  be  calculated  from  the  di¬ 
mensions  of  the  foils,  whereas  Rs  us  usually  unknown. 

We  can  also  define  the  effective  foil  impedance  [Q] 


which  contains  all  the  tab  geometry  information. 


2.4.  Foil  currents  and  potentials 

We  can  obtain  an  expression  for  the  foil  currents  by  using  (5) 
and  (6)  to  eliminate  the  stack  current  in  (10)  which  can  then  be 
integrated  to  give: 

lA(X)=^-Irc(2'c-Z'A)  (19) 

Ic(X)  =  -  IrA  (S'c  -  2J0  (20) 

Integrating  once  more  yields  explicit  expressions  for  the  foil 
potentials 

MX)  =  -LRa  j  IAdX  =  -rAcj>s(X)  +  IRfrArc(Sc  -  XA)  (21 ) 

MX)  =  -LRC  J  IcdX  =  rccf>s(X)  +  IRfrArc(Xc  -  SA)  (22) 

These  expressions  for  the  foil  potentials  are  important  for  two 
reasons,  1 )  they  allow  us  to  derive  explicit  expressions  for  the  cell 
voltage  and  impedance  under  load,  and  2)  they  are  used  to  deter¬ 
mine  the  fractional  currents  fj  for  complex  tab  configurations  by 
potential  matching: 


Rfoil  =  Rcell  “ 

=  Rf  (; rcr(Xc )  +  rAr(XA)  +  rArc[S(Xc)  -  X(X„)]  - 1) 

(27) 

foil  fraction  of  the  impedance  as 

r  .  =  Rfoil  =  j _ 2 _ 

f0"  Rceii  X2(rcr(Xc)  +  rAr(Xc)  +  rArc[E(Xc)-E(XA)}) 

We  need  to  be  a  little  careful  to  define  the  precise  meaning  of  Rs 
and  Rfoii.  Rs  only  corresponds  to  the  total  effective  stack  impedance 
if  the  stack  current  density,  is  is  uniform,  and  this  can  only  happen 
for  a  cell  with  continuous  current  collection.  As  soon  as  is  becomes 
non-uniform  the  effective  impedance  of  the  stack  becomes  slightly 
larger  than  Rs.  This  excess  is  actually  captured  in  the  definition  of 
Rfoii  in  Eq.  (27).  So  the  term  Rf0ii  is  a  little  bit  misleading.  Indeed  Rf0ii 
is  mostly  foil  impedance,  but  also  contains  a  small  contribution 
from  the  effective  stack  impedance  enhancement  caused  by  non- 
uniform  is.  An  accurate  name  might  be  Renhancement  over  continuous 
current  collection  ceil-  This  distinction  only  becomes  important  when 
we  disentangle  the  stack  and  foil  contributions  to  the  heat 
generation. 

The  total  heat  generation  rate  [watts]  in  the  cell  is 
1 

*7 cell  =  Qstack  +  Qfoil  =  l  J(wi2s(X)ps  +  I2a(X)Ra  +  /2(X)RC)dX 
0 


-  0a(^/<+ i)> 
0c(*k)  =  0c(W5 

2.5.  Cell  level  properties 


'  which  includes  stack,  anode  foil  and  cathode  foil  contributions. 

General  expressions  for  the  heat  are  rather  complex  and  will  not  be 
(24)  presented  here.  Expressions  for  specific  tab  configurations  will  be 
presented  below. 

In  Ref.  1]  we  described  a  method  for  modeling  spiral  wound 
cells  using  an  empirical  estimation  of  the  effective  foil  impedance. 
In  particular  we  approximated 


Defining  XA  as  the  location  of  any  anode  tab,  Xc  as  the  location  of 
any  cathode  tab  we  can  write  down  the  external  cell  voltage  as 


n  L(RA+RC) 

Rfoil  =  — 3g — 


(28) 
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where  g  is  the  so  called  geometry  factor  which  accounts  for  the  tab 
geometry.  Intuitive  arguments  were  employed  in  order  to  deduce 
ideal  values  of  g  that  were  valid  when  is  was  assumed  constant. 
More  accurate  values  of  g  were  determined  by  comparing  with  fully 
coupled  numerical  simulations.  We  can  combine  (28)  and  (26)  to 
obtain  an  explicit  expression  for  g 

g  =  -J - 1 - ry  (29) 

3  (rcr (Xc)  +  rAT(XA)  +  rArc[S(Xc)  -  2(X„)]  -  f) 


Table  2 

Three  special  states  during  a  constant  current  discharge. 


Time  into  discharge 

h(X) 

V0(X) 

Symbol  name 

Early 

Peaks  at  tabs 

Constant  in  X 

•early 

‘s 

Middle 

Constant  in  X 

Dips  at  tabs 

.•middle 

ls 

Late 

Dips  at  tabs 

Dips  at  tabs 

date 

ls 

inverted  and  is  peaked  between  the  tabs.  Our  definition  for  the  late 
stage  is 


which  is  valid  even  when  is  is  not  constant. 


date  _  n  .-middle  _  .-early 
ls  ~  ^ls  ls 


(33) 


2.6.  Time  evolution 

2.6.1.  Equation  of  motion 
By  re-arranging  Eq.  (2)  as 

psis(X,  t)  =  Vo(X,  t)  -  </v(X,  t)  +  IRfT(X)  (30) 

and  taking  the  time  derivative 


l 


o 


we  can  arrive  at  the  time  evolution  equation  for  the  stack  current 

l 

t(t)^  =  is (X,t)-X2  J  fa(X,S)is(S,t) ds  (31) 

0 

where  we  have  used 


Fig.  3  shows  what  the  stack  current  would  look  like  in  the  three 
states,  for  a  tab  configuration  with  two  cathode  tabs  at  X  =  {1,0}  and 
one  anode  tab  atX  =  {Vi}  (configuration  F4  below). 

The  theory  for  the  early  condition  is  straight  forward,  just  as¬ 
sume  a  constant  V0»  in  Eqs.  (15) — (17)  with  the  result  that  </v  =  V0» 
and  then  carry  this  result  through  to  the  foil  and  cell  properties.  The 
theory  for  the  middle  condition  contains  some  subtleties  and  is 
outlined  in  Appendix  A.  Since  the  cell  voltage  and  impedance  are 
both  linearly  related  to  is,  the  late  state  solution  for  these  cell 
properties  simply  follows  same  rule  as  used  in  Eq.  (33),  for  example 


plate 

Kcell 


2  R 


middle 

cell 


pearly 
(zell  • 


(34) 


Consideration  of  the  three  states  of  charge  are  only  important 
for  duty  cycles  that  use  the  cell  on  a  more  or  less  continuous  bases 
(constant  current  or  closely  spaced  pulses).  Readers  who  are  pri¬ 
marily  interested  in  pulsed  applications  for  which  the  cell  is 
allowed  to  relax  between  pulses,  will  only  need  to  consider  the 
early  stage  results  which  are  applicable  at  any  DOD  provided  V0  and 
ps  are  defined  at  that  DOD. 


9Vo 

at 


a Vq  dq 
dq  td t 


dVo 

dq 


*s(X,t) 


(32)  2.7.  Summary 


and  defined  the  time  constant 


9<7  |t 

Eq.  (31)  is  called  a  Barbashin  equation  [19,20]  with  stationary 
kernel  </>g-  What  these  equations  tell  us  is  that  if  the  system  starts 
out  with  a  Vq(X)  that  is  uniform  in  X,  and  is(X)  has  cusps  at  the  tabs; 
then  as  time  evolves  is(X)  gets  smoother  and  V0(X)  develops  peaks. 
For  constant  z  there  may  be  hope  of  finding  analytic  solutions  to 
(31 )  but  in  the  general  case  we  suspect  that  numerical  integration/ 
simulation  will  be  required. 


2.6.2.  Three  states  of  discharge 

Our  goal  in  this  work  is  to  describe  a  technique  for  comparing 
various  tab  configurations,  however  some  tab  configurations  will 
be  more  efficient  at  the  beginning  of  discharge  and  very  inefficient 
late  in  the  discharge  and  for  others  vice  versa.  So  we  need  an 
analytically  tractable  way  for  looking  at  the  solutions  and  cell 
properties  at  various  states  of  discharge  (DODs).  This  can  be  ach¬ 
ieved  by  considering  three  special  states  of  charge  summarized  in 
Table  2.  At  short  times  is(X)  has  cusps  at  the  tab  locations  and  V0  is 
constant,  somewhere  midway  through  a  discharge  is(X)  becomes 
roughly  constant  and  V0  had  downward  pointing  cusps  at  the  tab 
locations.  Finally  toward  the  end  of  a  discharge  is(X)  becomes 


Expressions  for  all  the  relevant  stack,  foil  and  cell  properties  are 
summarized  in  Table  3.  The  general  and  early  formulas  are  based  on 
hyperbolic  functions  in  X,  whereas  at  the  middle  state  the  formulas 
are  quadratic  in  X.  This  concludes  the  development  of  the  general 
theory  and  solutions.  In  the  remainder  of  the  article  we  will  discuss 
merits  of  specific  tab  configurations  and  a  number  of  examples 
showing  how  the  theory  can  be  applied. 
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Table  3 

Summary  of  all  formulas  for  general,  early  and  middle  state  conditions. 


Symbol 

General 

Early  t  =  0 

Middle 

Conditions 

None 

V0  =  const. 

is  =  const 

Stack  Potential  DE 

(V2  -A2)</>s  -  -X2V0+IR{Z" 

<=  with  const.  VQ 

v20s  =  IRf[Z"(X)  -  1] 

Particular  DE 

(V2  -  A2)</v  =  -X2V0(X) 

with  const.  VQ 

v2</v=  -iRf(i  -  A2r(x)) 

Particular  Solution 

MX) 

X2 /</>G(X,S)V0(S)dS 

/Kf[i  +  S(X)-±X2] 

Stack  potential  solution 

MX) 

M  -  iRf r 

V0  -  IRf r 

IRf[S(X)-XX2]  +  d 

Stack  current 

WO 

Vo-fa+iRf  r 

Ps 

I  r 

2LWA  1 

1 

2LW 

Anode  foil  current 

HX) 

- t-lrc( Zc-Za) 

7[_Vxr  -  rc(2'c  -  Z'A)\ 

I[2'a(X)  -X] 

Cathode  foil  current 

MX) 

_  irA(z'c  -  s'A) 

;[vxr  -  rA(E'c  -  S'A)] 

-/[X'C(X)-X] 

Anode  foil  potential 

MX) 

-rA<l>5+IRtrArc(Xc-2A) 

-IRfrA[XA(X)-^X2] 

Cathode  foil  potential 

MX) 

rc0s  +  IRf rArc(Xc  ~  2A) 

<= 

/Rfrc[Xc(X)-JX2]+d 

Cell  voltage 

vcen 

Eq.  (25) 

V0  —  TRcell 

Eq.  (A.4) 

Cell  impedance 

PceW 

Eq.  (26) 

Eq.  (26) 

Eq.  (A.5) 

3.  Results 

3.1  Internal  properties 

It  is  useful  to  look  at  all  the  internal  properties  for  a  simple  yet 
common  tab  configuration  with  one  anode  tab  at  X  =  0  and  one 
cathode  tab  at  X  =  1.  This  is  tab  configuration  F2  in  Table  5.  The 
details  of  the  cell  are  presented  in  Table  4.  We  purposely  chose 
asymmetric  foil  resistances  to  provide  understanding  of  this  effect. 
The  results  for  early,  middle  and  late  stages  of  the  discharge  are 
shown  in  Fig.  4.  The  following  features  should  be  noted: 

•  The  open  circuit  stack  potential  V0(X)  becomes  less  uniform  as 
the  discharge  proceeds. 

•  The  early  stack  current  is  peaked  at  the  tab  locations.  Because 
the  cathode  foil  is  higher  impedance  we  see  more  stack  current 
at  the  cathode  tab  than  at  the  anode  tab. 

•  For  the  early  stage,  near  the  cathode  tab  the  cathode  foil  current 
is  suppressed  and  the  anode  foil  current  is  enhanced.  The  sys¬ 
tem  prefers  to  run  as  little  current  as  possible  through  the 
higher  impedance  cathode  foil. 

•  The  middle  stage  foil  currents  are  linear,  which  is  just  a  conse¬ 
quence  of  the  constant  stack  current. 

•  At  the  late  stage  the  system  is  forced  to  run  more  foil  current 
near  the  cathode  tab  because  the  active  material  in  the  stack  has 
been  depleted  near  the  cathode  tab.  This  means  that  the  effec¬ 
tive  foil  impedance  of  the  cell  will  be  enhanced  at  the  late  stage. 

•  The  foil  potentials  look  smooth  and  fairly  similar,  but  this  is 
deceptive.  Small  changes  in  the  shape  and  value  of  the  foil  po¬ 
tentials  can  result  in  significant  changes  in  the  foil  currents  and 
stack  properties. 


Table  4 

Sample  cell  design  and  material  parameters. 


Property 

Units 

Value 

L 

[cm] 

80 

W 

[cm] 

6 

t 

[pm] 

10 

Pa,  Pc 

[pQ  cm] 

1.7,  3.1 

Pa,  Pc 

[mQ  cm-1] 

0.28,  0.52 

rA,  rC 

[no  units] 

1/3,  2/3 

Ps 

[Q  cm2] 

3.413 

X 

[no  units] 

3 

Vo 

[volt] 

4 

This  should  give  the  reader  some  insight  as  to  how  tabbed  cells 
prefer  to  distribute  current  and  why.  In  the  sections  that  follow  we 
will  focus  on  cell  impedance  for  systems  with  r&  =  rc  =  V2,  which  is 
what  most  commercial  cell  designs  strive  to  achieve.  This  restric¬ 
tion  will  greatly  simplify  the  formulas  presented  without  detracting 
from  their  practical  application.  This  also  means  that  for  a  specified 
tab  configuration,  all  results  are  equally  valid  if  the  anode  and 
cathode  tab  locations  are  swapped. 

3.2.  Cell  impedance  for  fixed  tab  configurations 

Table  5  shows  the  reduced  cell  impedance  (RCeii/^s)  and  geom¬ 
etry  factor  (g)  expressions  for  a  number  of  fixed  tab  configurations. 
For  each  configuration  we  show  a  drawing,  list  of  tab  locations, 
current  fractions  and  the  reduced  impedance  for  the  early  and 
middle  state  conditions.  The  late  stage  cell  impedance  can  easily  be 
calculated  from  Eq.  (34),  and  the  corresponding  geometry  factor 
should  be  calculated  from  using  (28).  However 

glate  =  2gmiddle  _  gearly  +  q ^  (35) 

is  valid  to  second  order  in  X  so  for  example  for  FI  g1310  =  1  -  A2/ 15 
which  points  us  towards  a  simple  method  of  defining  an  SOC 
dependence  for  g. 

Readers  can  simply  multiply  the  reduced  impedance  expres¬ 
sions  in  Table  5  by  either  Rf/A2  or  Rs  in  the  units  of  choice,  to  get  a 
real  cell  impedance.  For  most  configurations  the  fractional  currents 
were  easily  determined  by  uniqueness  (only  one  tab)  or  by  sym¬ 
metry.  For  configuration  F6  with  three  cathode  tabs,  the  potential 
matching  relations  (24)  were  used  to  derive  the  unsurprising 
answer  {1/4,  V2, 1/4}  for  the  cathode  tab  current  fractions.  Config¬ 
uration  F7  was  purposely  included  as  an  example  of  what  happens 
in  a  low  symmetry  configuration.  In  general  the  complexity  of  the 
expressions  increases  rapidly  as  the  symmetry  drops  or  the  number 
of  tabs  increases. 

For  the  early  state  geometry  factors  we  have  chosen  to  display  a 
second  order  Taylor  expansion  rather  than  the  full  expression.  For 
the  middle  state,  the  expressions  for  impedance  and  geometry 
factor  are  exact  (not  Taylor  expansions).  All  middle  state  imped¬ 
ances  are  independent  of  the  partial  foil  resistances.  Although  not 
evident  in  the  table  it  turns  out  that  the  same  is  true  for  the  early 
stage  impedance  for  tab  configurations  FI,  F3  and  F5.  This  means 
that  adjusting  the  relative  thicknesses  of  the  anode  cathode  foils 
will  yield  no  improvement  in  the  corresponding  cell  impedances 
for  these  three  tab  configurations. 
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Table  5 

Fixed  position  tab  configurations.  Drawings  use  red  for  anode  tabs  and  blue  for  cathode  tabs. 
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Name  Drawing 


FI 


F2 


F3 


F4 


F5 


F6 


F7 


Xj{A}{C]  fj{A}{C}  Early 

Rce\\ 

^stack 


{0}{o} 

Dili} 

{0}{i} 

{!}{!} 

{1H1} 

{!}{!} 

®{0,D 

aiHii} 

A  cosh  A 
sinh  A 


=  A  coth  A 


r 


A 

2  sinh  A 


(1  +  cosh  A) 


A  A 

2COth2 


^  H - - — ^fl+coshi) 

47  4  sinh 4  '  2  A 


A 

Tcoth- 
4  4 


(l  +cosh^ 


{1}{1  -  (C.l),  (C.l)} 


^(A  cosh  A  +  sinh  A) 
(cosh  A  -  1)  +  A  sinh  A 


Middle 

g  Keen  £ 

^stack 

;2  i 

1+T5  +  0(A4)  1+-A2  1 

1+^o  +  0(A4)  1+5a2  1 

4  +  ^  +  0(A4)  1+^A2  4 

4  +  240  +  0('A4')  1  +  12A2  4 

16 +  ^  +  0(A4)  1+^A2  16 

16  +  ^  +  0(A4)  1+^A2  16 

§+S+°(  a4)  1+^a2  ? 


There  are  a  number  of  interesting  patterns  in  Table  5: 

•  The  impedance  formulas  for  {FI,  F3,  F5}  make  a  nice  sequence. 
And  similarly  for  {F2,  F4,  F6}. 

•  To  lowest  order  the  geometry  factors  and  middle  stage  imped¬ 
ance  formulas  make  no  distinction  between  the  pairs  {FI,  F2}, 
{F3,  F4}  and  {F5,  F6} 

•  {FI,  F3,  F5}  are  all  strongly  dependent  on  A  (see  1/15  coefficient 
for  g)  for  the  early  and  late  stages.  Conversely  {F2,  F4,  F6}  are  all 
weakly  dependent  on  A  for  the  early  and  late  stages.  As  a  result 
we  would  expect  FI  to  have  a  steeper  voltage  curve  than  F2  but 
similar  average  voltage,  and  similarly  for  the  pairs  {F3,  F4}  and 
{F5,  F6}. 

There  are  a  number  of  interesting  graphical  comparisons  to  be 
made.  In  particular  FI,  F2  and  F3  all  have  two  tabs.  From  Figs.  5 
and  6  we  can  see  that  F3  is  vastly  superior  to  FI  and  F2  for  all 
states  of  charge  and  for  any  A.  Deciding  between  FI  and  F2  is  more 
subtle;  FI  is  superior  at  the  early  state  of  charge  but  inferior  to  F2  at 
the  late  stage.  Flence  FI  is  good  for  pulsed  applications,  but  F2  will 
have  a  flatter  voltage  curve  for  continuous  discharge  (energy)  ap¬ 
plications.  We  can  reach  the  same  conclusions  by  inspecting  the 
Taylor  expansions  for  g  in  Table  5  or  inspecting  the  coefficients  in 
Table  6.  From  Figs.  5  and  6  we  can  also  see  that  the  reduced  foil 
impedance  is  the  more  sensitive  value  for  comparing  various  tab 
configurations,  for  F4-F6  we  only  display  the  reduced  foil  imped¬ 
ance  in  Figs.  7  and  8.  The  general  trends  are  now  obvious,  there  is  a 
first  order  effect  for  impedance  {F5,  F6}  <  {F3,  F4}  <  {FI,  F2}  which 
is  controlled  by  the  first  term  in  the  Taylor  expansions  for  g.  This 
first  order  effect  overwhelms  any  dependence  on  state  of  charge  or 
A.  The  second  order  effects  tell  which  designs  will  have  the  flatter 
voltage  curves  {F2,  F4,  F6}  and  which  designs  are  better  for  pulsed 
applications  {FI,  F3,  F5}.  In  Fig.  9  we  show  the  geometry  factors  for 
F1-F6  and  also  the  foil  impedance  fractions  (Fig.  10).  The  foil 
impedance  and  geometry  factors  for  configurations  {F2,  F4,  F6}  look 


Fig.  4.  Internal  properties  for  early  (solid  line),  middle  (small  dashes)  and  late  (long 
dashes)  stages  of  discharge,  for  the  F2  tab  configuration. 
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0  2  4  6  8  10 

A  (No  Units) 


Table  6 

Taylor  expansion  coefficients  for  Eq.  (37). 


^0 

a2 

a4 

B2 

b4 

Max  X 

FI 

l 

3 

1 

“45 

2 

945 

0 

0 

4.1 

F2 

1 

1 

1 

1 

1 

10.3 

3 

720 

30,240 

12 

120 

F3 

F4 

1 

T2 

1 

1 

“720 

1 

1 

30,240 

1 

0 

1 

0 

1 

10.3 

26 

12 

11520 

1,935,360 

192 

7680 

F5 

1 

4S 

1 

11,520 

1,935,360 

0 

0 

26 

F6 

1 

48 

1 

184,320 

123,863,040 

1 

3072 

1 

491,520 

65 

F7 

5 

24 

19 

2880 

163 

483,840 

Undefined 

Undefined 

6.1 

g  ~  4.26  which  is  in  very  good  agreement  with  the  empirical  ad¬ 
justments  carried  out  against  the  numerical  simulation. 

In  Table  5  we  chose  to  express  the  reduced  cell  impedance  as 
Rcew/Rs  for  brevity  reasons,  however  Rs  is  not  a  directly  measurable 
quantity  and  it  is  difficult  to  calculate.  However  Rf  is  very  easy  to 
calculate  (see  Worked  examples  section  below).  Using  (12)  we  find 


^cell  _  J_  fRce\\\ 
Rf  ~a2\Rs) 


(36) 


Fig.  5.  Reduced  cell  impedance  for  the  FI,  F2  and  F3  tab  configurations  at  all  three 
stages  of  discharge. 


essentially  flat  and  to  a  good  approximation  we  can  ignore  the 
weak  A  dependence  for  these  cells. 

In  Ref.  [1  ]  we  used  Eq.  (28)  to  model  the  foil  impedance  of  a  cell 
and  compared  the  calculated  voltage  profile  under  load  with  a 
numerical  simulation,  using  a  large  enough  grid  to  give  very  ac¬ 
curate  results.  For  a  cell  design  with  the  F4  tab  configuration  and 
X  ~  8  we  found  that  adjusting  g  from  the  ideal  value  of  4  to  g  =  4.2 
gave  good  agreement  between  the  two  calculations.  For  foil  heating 
we  found  a  slightly  larger  value  g  =  4.3  gave  the  best  results.  From 
Table  5  we  see  that  for  this  design  g  =  4  +  A2/ 240  which  yields 


which  allows  one  to  express  the  cell  impedance  in  terms  of  Rf  and 
the  expressions  in  Table  5.  This  is  only  useful  of  we  already  know  A 
which  is  usually  not  the  case.  The  left  and  side  of  (36)  is  directly 
measurable  so  we  must  invert  (36)  in  order  to  determine  A.  This  is 
most  easily  achieved  with  a  Taylor  expansion 

%1  =  P+%l  =  F+^  +  (1-  2D0D)  (^a2  +^a4) + °(a6) 

(37) 

where  the  coefficients  for  all  tab  configurations  are  listed  in  Table  6. 
As  is  evident  from  (37)  by  simply  dropping  the  1/A2  term  at  the 
beginning  of  the  expansion  we  can  use  the  same  coefficients  to 
calculate  Rf0ii/Rf  if  desired.  In  Table  6  we  also  indicate  the  maximum 


0  2  4  6  8  10 

A  (No  Units) 


0  2  4  6  8  10 

A  (No  Units) 


Fig.  6.  Reduced  foil  impedance  for  the  FI,  F2  and  F3  tab  configurations  at  all  three  Fig.  7.  Reduced  foil  impedance  for  the  F3  and  F4  tab  configurations  at  all  three  stages 
stages  of  discharge.  of  discharge. 
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A  (No  Units) 


Fig.  8.  Reduced  foil  impedance  for  the  F5  and  F6  tab  configurations  at  all  three  stages 
of  discharge. 


value  of  A  for  which  the  expansion  (37)  is  accurate.  Max  A  is  defined 
as  the  point  where  the  sixth  order  term  in  (37)  becomes  of  order 
one.  Eq.  (37)  is  the  key  to  disentangling  the  foil  and  stack  contri¬ 
butions  to  the  total  cell  impedance.  All  we  need  is 

1.  A  measured  cell  impedance  (which  is  a  routine  measurement). 
This  could  be  an  AC,  DC  or  pulse  impedance. 

2.  Some  basic  information  on  the  foil  dimensions  to  calculate  Rf  =  L 

lPAl{WAtA)  +  Pcl(Wctc)] 

3.  The  tab  configuration  for  deciding  the  relevant  set  of  expansion 
coefficients  Aq,  A2  and  A4. 


20 
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_  10 
(/> 
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Is  10 

LL 
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g  20 
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15 
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0  2  4  6  8  10 

A  (No  Units) 

Fig.  9.  Geometry  factor  for  the  F1-F6  tab  configurations  at  all  three  stages  of 
discharge. 
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Once  A  is  determined  accurate  estimates  for  all  other  internal 
and  external  cell  properties  can  be  calculated.  See  worked  problem 
1  below. 

3.3.  Cell  heating  for  fixed  tab  configurations 

The  total  heat  generation  for  the  early  and  middle  stages  of 
discharge  can  be  simply  calculated  from  /2ftCeib  however  this  is  not 
true  of  the  late  stage  total  heat  or  for  the  breakdown  in  terms  of 
stack  and  foil  contributions.  Expressions  for  the  early  middle  and 
late  stage  foil  heating  rates  are  shown  in  Table  8  in  reduced  units  q/ 
I2Rs.  For  the  stack  only  the  early  stage  expressions  are  shown,  it 
turns  out  the  that  the  late  stage  makes  the  same  stack  heat  as  the 
early  stage  as  is  intuitively  evident  from  the  symmetry  in  Fig.  3.  The 
middle  stage  stack  heat  in  reduced  units  is  just  1. 

Before  discussing  the  results  and  trends  it  is  worth  under¬ 
standing  why  in  general  qf0ii  =£  I2Rm\-  This  is  not  because  of  ap¬ 
proximations.  As  discussed  above  in  Section  2.5  by  definition  Rf0[\  is 
mostly  foil  impedance,  but  also  contains  a  small  contribution  from 
the  effective  stack  impedance  enhancement  caused  by  non- 
uniform  is.  The  proper  interpretation  of  /2ftf0ii  is:  The  excess  heat 
generated  by  not  having  continuous  current  collection.  Fig.  12 
clearly  shows  that  for  the  FI  tab  configuration,  when  A  >  1  the 
heat  calculated  form  I2Rm\  is  slightly  overestimated  for  the  early 
stage  and  slightly  underestimated  for  the  late  stage.  The  same  is 
true  for  the  F3  and  F5  configurations.  For  the  F2  configuration  we 
see  that  this  distinction  is  tiny,  and  the  same  is  true  for  the  F4  and 
F6  configurations. 

Fig.  13  shows  the  heating  breakdown  for  all  three  states  of 
discharge  and  for  the  FI  and  F2  tab  configurations.  As  with  the  cell 
impedance,  we  see  that  for  the  heating  of  the  FI  configuration,  the 
dependence  on  the  stage  of  the  discharge  becomes  greatly  ampli¬ 
fied.  FI  is  not  good  for  full  discharge  applications,  but  good  for 
pulsed  if  it  has  time  to  relax  between  pulses.  The  stack  heating  may 
look  small  but  that  is  only  because  we  chose  to  normalize  by 
dividing  by  Rs.  If  we  had  normalized  with  Rf  we  would  see  the  stack 
heat  diverge  as  A  -►  0.  The  breakdown  for  the  other  four  high 
symmetry  tab  configurations  is  very  similar,  with  splitting  taking 
place  at  large  values  of  A. 

3.4.  Parametrized  tab  configurations 

In  this  section  we  will  describe  5  tab  configurations  for  which 
the  locations  of  some  or  all  of  the  tabs  are  determined  by  a  single 
parameter,  X.  Again  we  will  only  present  expressions  for  the  case 
rA  =  rc  =  V2.  It  is  tempting  to  minimize  these  impedance  expressions 
in  order  to  find  optimal  tab  locations,  but  this  is  only  useful  in  low 
symmetry  configurations  not  presented  here.  For  the  configura¬ 
tions  in  Table  7  the  optimal  values  X  is  usually  evident  by  inspec¬ 
tion.  The  reduced  foil  impedances  for  all  parametrized 
configurations  are  shown  in  Fig.  11  where  we  see  that  the  tab  lo¬ 
cations  used  in  Table  5  are  all  optimal.  The  parametrized  expres¬ 
sions  are  probably  most  useful  when  limitations  of  cell  design  or 
manufacturing  equipment  forces  one  to  deviate  from  the  optimal 
tab  locations.  The  expressions  in  Table  7  can  then  be  used  to 
calculate  the  relative  or  absolute  cost  in  terms  of  increased  cell 
impedance  caused  by  deviations  from  the  optimal  locations.  These 
expressions  can  also  be  used  to  disentangle  the  effective  foil 
impedance  from  stack  ASI  for  many  low  symmetry  tab 
configurations. 

4.  Worked  examples 

In  this  section  we  will  present  some  examples  of  potential  real 
world  problems  that  can  be  easily  solved  using  the  results  derived 
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Table  7 

Parametrized  tab  configurations. 


Name  Drawing 


PI 


Xj{A}{C) 

fmc] 

Early 

^stack 

Middle  ^cel1 

^stack 

{X}{X} 

{!}{!} 

X  cosh(AX)cosh(AX) 
sinh  X 

1  +1A2(1  -3XX) 

TOW 

{!}{!} 

C.2 

1  +1A2(1  -3XX) 

{X,X}{0,1} 

C.3 

{X,X}{X,X} 

C.4 

^+X2[X(X-\)+^-H\X-\)\ 

{i  +  X}{0,l} 

{1}{C.5,  1  -  C.5} 

C.6 

1+^A2(5-12XX) 

above.  Each  example  will  be  presented  in  the  form  of  a  question, 
supplemented  by  some  required  details  and  then  the  answer  or 
solution. 


3.  Calculate  Rcei\lRf  =  0.785  [no  units] 

4.  Using  the  coefficients  for  the  F2  configuration  in  Table  6, 
DOD  =  0  and  Eq.  (37)  we  build  the  equation  we  need  to  solve 


4.1.  I  have  measured  the  impedance  of  my  cell.  How  do  I  subtract 
the  effective  foil  impedance  in  order  to  obtain  the  correct  stack  ASI? 


0785  *JE+5-H0J2  +  30540 


-A 


(38) 


Details :  Anode  foils  are  83  cm  x  6.2  cm  x  10  pm,  cathode  foils 
are  70  cm  x  6.0  cm  x  19  pm.  There  is  a  cathode  tab  at  the  core,  and 
an  anode  on  the  outside  of  the  jelly  roll.  The  tabs  are  right  next  to 
where  the  active  coatings  start.  The  measured  cell  impedance  was 
30  [mQ] 

Solution :  This  is  a  very  common  problem  which  to  our  knowl¬ 
edge  has  not  been  solved  in  literature  before.  We  will  use  the  well 
known  material  properties:  1.7  [pQ  cm]  and  3.1  [pQ  cm]  for  the  Cu 
and  A1  resistivities.  Based  on  the  description,  this  cell  has  the  F2  tab 
configuration  with  the  anode  and  cathode  tabs  reversed,  which  has 
no  effect.  The  steps  are: 

1.  Calculate  Ra  =  paK^a^a )  =  1.7  x  10-6  Q  cm/(6.2  cm  x 
0.001  cm)  =  0.274  [mQ  cm-1]  and  Rc  =  3.1  x  10“6  Q  cm/ 
(6.0  cm  x  0.0019  cm)  =  0.272  [mQ  cm-1  ].  We  are  in  luck  RA  ~  Rc 
so  we  need  not  worry  about  the  effects  of  the  fractional  foil 
impedances  rA,  rc. 

2.  Use  I  =  70  cm  (not  83  cm)  to  calculate  Rf  =  L(RA  +  Rc)  =  38.2 
[mQ] 


Table  8 

Analytic  expressions  for  the  foil  and  stack  heating  rates. 

Name  ^ _  wk_ 

Early  Middle  Late  Early 


FI 

F2 

F3 

F4 

F5 

F6 


K 


isinh  2A-AN 
cosh  2A-1 ) 


A2 

~4 


+ 


A 

4 


( 


sinh  A — A  \ 
cosh  A-l/ 


A 

4 


( 


sinh  A-A\ 
cosh  A  1 ) 


|A2 

J*2 


4  +  <  +  C'y-4cothA  A(if^) 


4  +  ^  +  <fi‘y-2 


A_  A  (  sinh  A+A^\ 

1  4  ^cosh  A  1 ) 

a  |  A2  |  nearly  ->  A  sinh  A  A  ( sinh  A+A^\ 

4  +  T  +  qfoil  -^cosh  4  i  cosh  A-l ) 


A2  ,  A 
T6  +  T6 


f  2sinhA-A\ 
cosh  A-l  J 


4  +  ft  +  C'y- 


A  sinh  A 
cosh A-l 


A  ( 2sinhA+A\ 
16  cosh  A-l  J 


±  ( 2sinhA-A\ 
16  cosh  A-l  J 


1  A2 


4  +  h  +  <Im 


early 


A  sinh  A 
cosh  A-l 


( 2  sinh  A+a\ 
16  cosh  A-l  J 


A2  ,  A 
64  +  64 


4sinhA-A\ 
cosh  5- 1  J 


1  ;2 
48X 


4  1  A2  ,  yearly 
4  +  48  +  4foil 


1  A  sinh  A 

2  cosh^-1 


_A_  / 4sinhA+A\ 
64  ycosh^-1  J 


5.  Using  a  numerical  solver  we  find  A2  ~  2.200  [no  units].  This 
result  is  much  less  than  max  A  =  10.3  so  we  can  be  confident  that 
the  expansion  is  accurate. 


0  0.2  0.4  0.6  0.8  1 

x/L  (No  Units) 

Fig.  11.  Reduced  foil  impedance  for  the  P1-P5  parametrized  tab  configurations  at  all 
three  stages  of  discharge. 
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A  (No  Units) 
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Fig.  12.  Foil  heating  from  the  exact  calculation  and  from  I2Rm\ 


6.  Use  Eq.  (12)  to  find  Rs  =  Rf/A2  =  17.4  [mfi] 

7.  Use  Eq.  (12)  to  find  the  stack  ASI  ps  =  2 LWRS  = 
2  x  70  cm  x  6  cm  x  0.0174[Q]  =  14.6[Q  cm2]. 


A  (No  Units) 


Fig.  13.  Foil,  stack  and  total  heat  generation  rates  for  the  FI  and  F2  tab  configurations. 


Anode  1cm 


Cathode 


2cm 


Fig.  14.  Electrode  design  with  bare  foil  regions  near  the  tabs,  as  described  in  worked 
example  #5. 


If  we  had  simply  assumed  that  the  foil  impedance  was  negligible 
and  used  ps  =  2LWRce\\  =  25.2  [0  cm2]  our  result  would  be  over¬ 
estimated  by  about  70% !  In  general  once  you  know  A  everything 
becomes  “unlocked”  and  you  can  calculate  all  the  internal  and  cell 
level  properties  describe  above.  The  foil  impedance  is  relatively 
insensitive  to  frequency  so  if  your  measured  cell  impedance  was 
from  a  DC  pulse  power  test  then  the  derived  value  for  ps  will  be  the 
DC  stack  ASI,  and  correspondingly  for  AC. 

Usually  one  can  get  good  results  with  lower  order  versions  of  Eq. 
(38).  For  example  if  we  ignore  the  last  term  and  solve 


Jb 2 

720 


(39) 


the  result  is  A2  =  2.199  [no  units],  i.e.  the  difference  is  in  the  third 
decimal  place.  In  fact  even 


0.785  =  T  +  I 

A2  3 


(40) 


gives  a  reasonable  result  A2  ~  2.214  [no  units].  The  bottom  line  is 
that  for  most  problems  you  don’t  need  the  higher  order  terms. 


4.2.  My  cell  has  an  F3  tab  configuration  but  I  am  out  of  luck, 

Ra  =£  Rcso  I  cannot  use  the  tidy  formulas  in  Table  5 

Solution:  On  the  contrary,  you  are  in  luck,  the  formulas  for  the  FI, 
F3  and  F5  tab  configurations  just  happen  to  be  independent  of  the 
fractional  foil  impedances  rA,  rc  so  Table  5  is  still  valid. 


4.3.  My  cell  has  an  F4  tab  configuration  and  now  I  am  really  out  of 
luck,  Ra  =£  Rcso  I  cannot  use  the  tidy  formulas  in  Table  5.  How  do  I 
calculate  the  stack  ASI? 

Details:  Anode  foils  are  83  cm  x  6.2  cm  x  15  pm,  cathode  foils 
are  70  cm  x  6.0  cm  x  15  pm.  The  measured  cell  impedance  was 
13  mQ 

Solution:  Steps  1—3  are  the  same  as  for  problem  #1  above. 

1.  Ra  =  0.183  mQ  cm-1  and  Rc  =  0.344  mQ  cm-1 

2.  rA  =  0.347,  rc  =  0.653 

3.  A r2  =  0.0235 

4.  Rceii/Rf  =  0.352 

5.  Using  the  coefficients  for  the  F4  configuration  in  Table  6  and  Eq. 
(B.l)  we  build  the  equation  we  need  to  solve 
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1935360 
1 


002351  -T^+TSSO1' 


6.  Using  a  numerical  solver  we  find  A1 2 3 4 5  ~  3.712. 

If  we  had  ignored  the  asymmetry  (i.e.  assume  A r  =  0)  we  would 
have  obtained  A2  ~  3.718!  In  general  it  is  not  easy  to  contrive  a 
realistic  cell  design  for  which  Ar  is  actually  important. 

4.4.  How  do  I  account  for  foil  impedance  in  my  cell  model? 

Solution: 

1.  Work  out  A  using  ps  calculated  from  the  stack  model.  ps  can  be 
calculated  analytically  as  described  in  the  appendix  of  Ref.  [1  ],  or 
by  doing  a  finite  difference  estimation  ps  ~  A0s/Ais. 

2.  Choose  an  expression  for  g  for  your  tab  configuration  from 
Table  5  and  use  A  to  calculate  a  numerical  value  for  g.  For  full 
discharges  you  can  make  g  DOD  dependent  using  Eq.  (35)  and 
follow  (37). 

3.  Subtract  IRfl(3g)  from  the  calculated  cell  voltage  in  your  model. 

4.  Add  I2Rfl(3g)  to  the  generated  heat  in  your  model.  Use  the  ex¬ 
pressions  in  Table  8  if  you  need  a  breakdown  in  terms  of  stack 
and  foil  contributions  to  the  total  heat. 


Fig.  10.  Foil  fraction  of  the  cell  impedance  for  the  F1-F6  tab  configurations  at  all  three 
stages  of  discharge. 


5.  Plot  IA{X),  Jc(X),  un  0c(X),  0s(X)f  2s(X)  using  (19),  (20),  (21), 
(22),  (17),  (2)  so  you  can  see  what  is  going  on  as  the  simula¬ 
tion  proceeds. 

This  is  not  the  only  way  to  use  these  techniques  to  enhance  a  cell 
model,  but  it  is  the  simplest  method. 

4.5.  My  cell  has  bare  foil  segments  near  the  tabs ,  your  derived 
expressions  in  Tables  5  and  7  don’t  account  for  that,  so  guess  I  am 
out  of  luck? 

Nope.  There  is  no  stack  current  and  therefore  no  foil  current 
divergence  in  the  bare  foil  regions  so  you  can  just  add  the  imped¬ 
ance  of  these  segments  to  your  hardware  impedance. 

Details :  Anode  foils  are  180  cm  x  6.2  cm  x  10  pm,  cathode  foils 
are  170  cm  x  6.0  cm  x  19  pm.  F4  tab  configuration  with  anode  and 
cathode  tabs  swapped.  1  cm  bare  foil  on  each  side  of  the  anode  tabs, 
2  cm  bare  foil  on  each  side  of  cathode  tab.  See  Fig.  14. 

Solution :  Any  bare  foil  between  tabs  and  the  end  of  foil  is 
irrelevant,  there  are  no  foil  currents  flowing  there. 

1.  Ra  =  0.274  mQ  cm-1  and  Rc  =  0.272  mQ  cm-1. 

2.  For  the  F4  tab  configuration  we  need  to  worry  about  current 
splitting  when  calculating  the  effective  impedance  of  the  bare 
foil  regions.  The  foil  currents  will  be  7/2  in  all  bare  regions  for 
this  cell. 

3.  RH  =  ^header  +  Rweids  +  Ra'M  1  cm  +  1  cm)  +  RCV i{2  cm  +  2  cm) 

4.  If  you  are  modeling  subtract  /Rh  from  the  calculated  cell  voltage 
and  add  /2Rh  to  the  generated  heat  in  your  model. 

5.  If  you  are  calculating  the  stack  ASI  as  in  problem  1,  then  subtract 
Rh  from  your  measured  cell  impedance  before  beginning  the 
steps  in  problem  1. 


This  will  be  very  accurate  for  tab  configurations  like  {F2,  F4,  F6} 
and  usually  “good  enough”  for  other  tab  configurations.  See  Ref.  [1  ] 
for  some  examples. 

4.7.  I  want  to  understand  the  impedance  of  a  cell  with  a  large 
number  of  tabs,  do  I  need  to  resort  to  numerical  methods  in  this  case 
because  the  algebra  will  get  too  complicated ? 

Solution:  Not  necessarily.  For  example  Lee  et  al.  [21]  did 
large  scale  numerical  simulations  for  cells  with  2,  5,  10  and  20 
tabs  per  electrode  at  locations  XA  =  Xc  =  {1  /TV,  3/N,  5/N...N  -  1/ 
N }  where  N  is  the  total  number  of  tabs.  They  called  these 
configurations  T2,  T5,  T10  and  T20  respectively.  Their  T2 
configuration  is  the  same  as  the  F5  configuration  in  Table  5. 
These  are  high  symmetry  configurations  and  as  a  result  we  can 
derive  an  expression  for  the  reduced  the  cell  impedance  for 
general  N 


pearly 

Kcell 


A 

=  N  C°th  N 


ggj l!  =  w2 

Rs  3\n)  ■ 


Nothing  overly  complicated  about  those  expressions !  Lee  et  al. 
calculated  voltages  curves  which  show,  what  looks  like  constant 
(ohmic)  shifting  with  N.  They  could  have  got  almost  exactly  the 
same  result  with  a  simple  lumped  foil  model  as  described  in 
problem  4  above. 


4.6.  I  don’t  know  how  to  use  numerical  solvers  and  I  just  want  a  4.8.  I  want  to  follow  the  steps  in  problem  1,  but  my  tab 

rough  estimate  of  the  foil  impedance  in  my  cell  configuration  is  not  in  any  of  the  tables,  now  what? 


Solution: 

•  Just  use  the  middle  state  geometry  factors  in  Table  5  and  then 
plug  g  into  Rfoii  =  Rf/(3gJ. 


Solution:  If  your  tab  configuration  is  not  in  the  table  then  is  it 
likely  to  be  a  low  symmetry  configuration  with  multiple  tabs  per 
electrode,  which  means  that  the  fractional  currents  cannot  be 
determined  by  symmetry  or  uniqueness.  Flere  are  the  steps: 
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1.  Brace  yourself! 

2.  Convert  your  tab  locations  to  dimensionless  units  using  the 
active  length  L 

3.  Use  Eq.  (18)  to  calculate  expressions  for  T (Xj)  for  all  tab 
locations. 

4.  Calculate  expressions  for  Sc{Xj)  -  SA(Xj)  using  Eqs.  (8)  and  (7) 
with  the  delta  functions  replaced  by  ramp  functions  (see 
Table  1 ).  Do  this  for  all  tab  locations. 

5.  Calculate  fA{Xj)  for  all  anode  tab  locations.  Do  the  same  for 
cathode. 

6.  Now  solve  the  potential  matching  Eqs.  (23)  and  (24)  and  extract 
expressions  for  the  fractional  currents. 

7.  Choose  any  anode  tab  location  XA  and  any  cathode  tab  location 
Xc.  If  the  potential  matching  was  done  correctly  it  should  make 
no  difference  which  tabs  you  choose. 

8.  Calculate  rcT(Xc)  +  rAr(XA)  +  rArc[S(Xc)  -  X(XA)]  which  will  be 
some  messy  function  of  A. 

9.  Use  Eq.  (26)  to  solve  for  A  as  per  problem  1. 

We  suggest  you  either  use  a  computer  algebra  system  or  email 

the  corresponding  author  and  he  will  do  it  for  you. 


5.  Conclusions 

The  analytic  theory  for  ID  foil  currents  and  potentials  in  spiral 
wound  cells  with  arbitrary  tab  configurations  has  been  derived.  The 
theory  has  been  extended  to  approximately  account  for  accumu¬ 
lation  of  the  gradients  in  the  state  of  charge  along  the  length  of  the 
foils.  Analytic  expressions  for  cell  level  properties  were  derived  for 
specific  tab  configurations.  The  quality  of  many  tab  configurations 
was  compared.  We  found  that  the  common  tab  configurations  can 
be  categorized  by  their  dependence  on  the  dimensionless  foil  to 
stack  impedance  ratio,  A.  We  learned  that  for  tab  configurations 
with  opposing  anode  and  cathode  tabs  the  cell  impedance  is 
strongly  affected  by  A  and  state  of  charge  and  conversely  when  the 
anode  and  cathode  tabs  are  well  separated  the  A  dependence  is 
relatively  weak.  It  also  became  clear  after  working  a  realistic  nu¬ 
merical  example  that  asymmetries  between  anode  and  cathode  foil 
conductivities  will  have  a  very  minor  effect  on  the  overall  cell 
impedance,  but  do  have  some  effect  on  the  current  distributions.  A 
straight  forward  method  is  described  for  disentangling  the  foil  and 
stack  components  of  a  measured  cell  impedance.  The  method  re¬ 
quires  one  to  effectively  invert  a  non  linear  function  of  A  thereby 
showing  why  this  useful  process  has  never  been  obvious.  The  an¬ 
alytic  theory  shows  us  how  to  directly  calculate  the  foil  geometry 
factor  introduced  in  a  previous  article  [  1  ]  and  we  find  the  results 
are  in  excellent  agreement  with  an  empirically  adjusted  value  used 
in  numerical  simulations.  From  the  analytic  results  it  came  to  light 
that  some  extra  heat  is  generated  in  the  stack  resulting  from  the 
non-uniform  stack  current.  It  is  a  matter  of  semantics  and  taste  as 
to  whether  or  not  one  includes  this  as  part  of  the  foil  heating.  We 
also  describe  a  simple  method  for  adding  foil  effects  to  a  stack¬ 
centric  cell  model.  It  should  be  clear  from  this  article  that  analytic 
results  provide  a  lot  of  insight  and  ease  of  use  that  numerical  so¬ 
lutions  cannot  deliver.  Readers  should  be  aware  this  theory  does 
not  account  for 

1.  Complex  voltage  curves.  This  will  only  affect  the  estimates  of  the 
middle  and  late  stage  cell  properties. 

2.  Non-uniformities  in  the  stack  resistivity  ps.  Again  this  is  only 
important  for  the  estimates  of  the  middle  and  late  stage  cell 
properties. 

3.  Non-ohmic  stack.  Yet  again,  this  is  only  important  for  the  esti¬ 
mates  of  the  middle  and  late  stage  cell  properties 


4.  2D  foil  effects  [6  .  These  would  be  important  in  a  very  tall  cell 
(large  W)  with  edge  joined  tabs. 

With  these  caveats  in  mind  we  would  suggest  that  cell  models 
can  still  get  very  accurate  results  using  the  solutions  derived  in  this 
article,  and  not  have  to  always  resort  to  numerical  solutions  for  the 
foil  portion  of  the  model. 

The  most  important  item  on  the  agenda  for  future  work  is  to 
extend  this  analytic  theory  to  handle  2D  systems  which  are  rele¬ 
vant  for  flat  pouch  cells.  Preliminary  work  on  this  problem  indicates 
that  it  is  also  analytically  tractable  but  very  challenging. 

List  of  symbols 


For  each  symbol  we  provide  example  units  in  order  to  define  the 
dimensions. 

Ai  Taylor  expansion  coefficients  for  the  cell  impedance  [no 

units] 

DOD  depth  of  discharge,  ranges  from  0  to  l[no  units] 

/foil  fraction  of  the  total  cell  impedance  caused  by  the  effective 

foil  impedance  [no  units] 

fj  fraction  of  the  total  external  current  flowing  into  tab  j  [no 

units] 

g  foil  impedance  geometry  factor  [no  units] 

H(X)  Heaviside  step  function  [no  units] 

H^X)  =XH(X),  ramp  function  [no  units] 

is(x)  current  density  across  one  layer  of  the  stack  [amp  cm-2] 

I  total  external  current  applied  to  the  cell  [amp] 

1/ 1,  Ic  anode  and  cathode  foil  currents  [amp] 

L  active  length  of  the  foils  [cm] 

nA,  nc  number  of  anode  and  cathode  tabs  [no  units] 

N  =nA  +  no  total  number  of  tabs. 

q  charge  [Ah] 

qCeii  total  heat  generation  rate  in  the  cell  [watt] 
qfoii  heat  generation  rate  from  the  foil  currents  [watt] 

Qstack  heat  generation  rate  from  the  stack  current  [watt] 

RA,  rc  =pA/WAtA ,  pc/Wctc,  linear  anode  and  cathode  foil 

resistivities  [Q  cm-1] 

rA,  rc  fractional  foil  impedances,  rA  +  rc  =  1  [no  units] 

Ar  deviation  of  rA  and  rc  from  Yi 

Rceii  total  cell  impedance  [Q] 

Rfoii  effective  impedance  of  the  foils  [Q] 

Rf  =L(Ra  +  Rc),  characteristic  impedance  of  the  foils  [Q] 

Rs  =ps/2LW,  total  resistance  of  the  stack  [Q] 

S  convolution  integration  variable  [no  units] 

t  Time  [s] 

tA,  tc  thicknesses  of  anode  and  cathode  foils  [pm] 

V0{x)  open  circuit  voltage  of  the  stack  [volt] 

Vceii  cell  voltage,  possibly  under  load  [volt] 

Vqcv  open  circuit  voltage  of  the  cell  [volt] 

WA,  Wc  widths  of  anode  and  cathode  foils  [cm] 

W  mm{WA,  Wc]  active  width  [cm] 

x  distance  along  the  length  of  the  foils  [cm] 

y  distance  across  the  stack  [pm] 

z  distance  across  the  width  of  the  foils,  or  the  height  of  the 

cell  [cm] 

X  =x/L,  dimensionless  distance  along  the  length  of  the  foils 

[no  units] 

X  =1  -  X,  mirror  of  distance  along  the  foils  [no  units] 

XA,  Xc  locations  of  any  anode  or  cathode  tab. 

V*  dimensioned  gradient  operator  a/6x  [cm-1] 

Vx  dimensionless  gradient  operator  9/aX  [no  units] 
a1 2 3  foil  to  stack  impedance  ratio  [cm-2] 
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<5(x)  Dirac  delta  function. 

(/>Af  (pc  anode  and  cathode  foil  potentials  [volt] 

0gPC  Xj)  Greens  function  for  unit  flux  at  X  =  Xj  [no  units] 

0s(x)  potential  difference  across  the  stack  [volt] 

0v(x)  open  circuit  portion  of  the  stack  potential  solution  [volt] 
T(X)  weighted  sum  of  Greens  functions  [no  units] 

A2  = o /2L2,  dimensionless  foil  to  stack  impedance  ratio  [no 

units] 

Pa,  Pc  anode  and  cathode  foil  material  resistivities  [Q  cm] 
ps  area  specific  stack  impedance  (ASI)  [0  cm2] 

t  time  constant  in  the  Barbashin  equation  of  motion  [s] 

Superscripts: 

early  early  stage  of  a  constant  current  discharge 
middle  middle  stage  of  a  constant  current  discharge 
late  late  stage  of  a  constant  current  discharge 

Subscripts: 

A  anode 

C  cathode 

j  tab  index 

s  stack 


Expressions  for  foil  currents  and  potentials  are  obtained  by 
directly  integrating  (5)  and  (6)  etc.  and  are  shown  in  Table  3.  We 
can  integrate  these  to  get 


Vcell  =  <£c(*c)  -  M*a) 

=  IRf  [rc2c(Xc)  +  rASA(XA)  -  X-fx\  - gxg]  +  d  (A.4) 


^cell  — 


Vocv  -  V, 


cell 


=  Rt 


i+E^2 


g  -  rc2c(Xc)  -  rAZA(XA)  +  X-fx\  +  Xj-X2c 

(A.5) 


Appendix  B.  Asymmetric  foil  resistance 

When  Ra  =£  Rc  then  we  cannot  assume  =  rc  =  Vi  as  was  done 
in  Table  5.  In  this  case  the  resulting  formulas  become  more 
complicated.  However  we  can  easily  extend  Eq.  (37)  to  handle  the 
asymmetric  case  by  setting  rA  =  Vi  +  Ar  and  rc  =  Vi  -  A r 


Appendix  A.  Middle  state  theory 

When  the  stack  current  is  constant:  is  =  I/2LW  =  constant  Eq. 
(30)  becomes 

V0(X)  =  4>S(X)+IRS.  (A.l) 

In  addition  Eq.  (10)  becomes 

Vjfo  =  «f[S"  -  1]  (A.2) 


^  ^  +  A0  +  (1  -  2DOD)  (a2X2  +  A4A4) 

+  Ar2(l  -  2DOD)  (b2X2  +  B4A4)  +o(A6)  (B.l) 

This  works  for  all  the  configurations  in  Table  5  except  F7. 

Appendix  C.  Long  formulas 


which  can  be  directly  solved  by  integration 


=  IRf 


2-J*2 


+  d. 


(A.3) 


The  integration  constant  d  is  related  to  stack  OCV.  As  mentioned 
above  V0  represents  the  fully  relaxed  open  circuit  voltage  across  the 
stack  long  after  current  is  turned  off.  If  we  take  the  average  this 
value  we  get  an  approximation  for  the  fully  relaxed  OCV  of  the  cell. 
The  average  can  be  carried  out  using  (A.l)  and  (A.3) 


^iocv  — 


1 

/  Vo(X) 


=  IRf 


+  d  +  IRS 


which  can  be  solved  for  d 


d  —  Vnry  —  IRf 


i  +  E  nfitf-l 


LAZ  ^  -  J  6J 

and  the  final  expression  for  the  stack  potential  is 


=  IRf 


fd  = 


cosh  A  -  1 


cosh  A  -  1  +2sinh  ^ 


A 


mn2 


16  1 


A 


+  si„h(J|2X-l|)]+rea5Fj 
x  (2  (cosh  X  +  l)(cosh(AX)  +  cosh  (AX)  +  1) 
+  [cosh(AX)  +cosh(AX)]2| 


(Cl) 


yearly  ^  _ |2X  -  !|  +  _I_]_I[cosh(AX)  +  COSh(AX)] 
-Asinh(A|2X-l|) 

l2  1  — 

i?ear*y  =  '6(1  |2X  l|)-1f.[2  sinh(AX)  +  2  sinh(AX) 


(C.2) 


(C.3) 


R6arly  =  4sinhA(c0Sh(AX)  +  cosh(AX))2  - ^sinh(A|2X -  1|) 


(C.4) 


/ci 


1  (cosh  A  -  1)  [2  -  cosh(AX)  +  cosh(AX)]  +  sinh  A[2AX  +  sinh(AX)  -  sinh(AX)] 

2  2(cosh  A  -  1)  +  A  sinh  A 


(C.5) 
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yearly  =  ^j  p AX- sinh(AX)  +  sinh(AX)]/cl  - sinh(AX)  - sinh(AX) } 

+ -  -  /  ^  j  { (cosh  A + 1 )  [cosh(  AX)  +  cosh  (AX)  ]  +  (cosh  A  - 1 ) 

x  [2  +  cosh(AX)-cosh(AX)]/cl+2 
+2cosh(AX)cosh(AX)} 

(C.6) 
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